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Abstract 

Equilibrium statistics of Hamiltonian systems is correctly described by the micro- 
canonical ensemble. Classically this is the manifold of all points in the A^— body 
phase space with the given total energy. Due to Boltzmann-Planck's principle, 
= tr{6{E — H)), its geometrical size is related to the entropy S{E, N,V, ■ ■ ■). 
This definition does not invoke any information theory, no thermodynamic limit, 
no extensivity, and no homogeneity assumption. Therefore, it describes the equi- 
librium statistics of extensive as well of non-extensive systems. Due to this fact 
it is the fundamental definition of any classical equilibrium statistics. It addresses 
nuclei and astrophysical objects as well. S{E, N,V, ■ ■ ■) is multiply differentiable 
everywhere, even at phase-transitions. All kind of phase transitions can be distin- 
guished sharply and uniquely for even small systems. What is even more important, 
in contrast to the canonical theory, also the region of phase-space which corresponds 
to phase-separation is accessible, where the most interesting phenomena occur. No 
deformed q-entropy is needed for equilibrium. Boltzmann-Planck is the only appro- 
priate statistics independent of whether the system is small or large, whether the 
system is ruled by short or long range forces. 
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1 Introduction 

Classical Tliermodynamics was originally designed to describe steam engines. 
It got its theoretical foundation in classical thermostatistics. Now Thermody- 
namics and the theory of phase transitions of homogeneous and large systems 
are some of the oldest and best established theories in physics. It may look 
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surprising to add anything new to it. Let me recapitulate what was told us 
since ~ 150 years: 

• Thermodynamics addresses large homogeneous systems at equilibrium (in 
the thermodynamic limit N oo\N/v=p,homogeneous) ■ 

• Phase transitions are the positive zeros of the grand-canonical partition 
sum Z(T,fj,, V) as function of e^'^ (Yang-Lee-singularities). Of course these 
singularities indicate the break-down of the (grand-) canonical ensemble. 

• Micro and canonical ensembles are equivalent. 

• Thermodynamics works with intensive variables T, P, /x. 

• Unique Legendre mapping T ^ E. 

• Heat only flows from hot to cold (Clausius). 

• Second Law only in infinite systems when the Poincare recurrence time 
becomes infinite (much larger than the age of the universe (Boltzmann)). 

Under these constraint only a tiny part of the real world of equilibrium systems 
can be treated. The ubiquitous non-homogeneous systems: nuclei, clusters, 
polymers, soft matter (biological) systems, but also the largest, astrophysical 
systems are not covered. Even normal systems with short-range coupling at 
phase separations are inhomogeneous and can only be treated within conven- 
tional homogeneous thermodynamics (e.g. in van-der-Waals theory) by bridg- 
ing the unstable region of negative compressibility by a Maxwell construction. 
Thus even the original goal, for which Thermodynamics was invented some 
150 years ago, the description of steam engines is only artificially solved. There 
is no (grand-) canonical ensemble of phase separated and, consequently, inho- 
mogeneous, configurations. This has a deep reason as I will discuss below: 
here the systems have a negative heat capacity c (resp. susceptibility). This, 
however, is impossible in the (grand-) canonical theory where c oc {SEy 

As was recently remarked by Pitowsky(0): "There is a schizophrenic attitude 
in the foundations of statistical mechanics. While Boltzmann's view has been 
promoted as conceptional superior to that of Gibbs, it is the canonical and 
not the microcanonical probability distribution that is extensively used in cal- 
culations. The claim being that the mathematical simplifications aside, all the 
ensembles are thermo dynamically equivalent because as long we are dealing in 
systems containing a large number of molecules.". Here, I will show that this is 
not only schizophrenic but is merely wrong in the most interesting situations. 
And that not only for small systems but also for the really large ones. 



2 Boltzmann-Planck's principle 

The Microcanonical ensemble is the ensemble (manifold) of all possible points 
in the 6N dimensional phase space at the prescribed sharp energy E: 
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W{E, N, V) = eotr6{E - H^) = e, j j^J^-^5{E - H^). 

Thermodynamics addresses the whole ensemble. It is ruled by the topology of 
the geometrical size W{E, N,- ■ ■), as expressed on Boltzmann's epitaph: 




which is the most fundamental definition of the entropy S. Entropy and with it 
micro-canonical thermodynamics has a pure mechanical, geometrical founda- 
tion. No information theoretical formulation is needed. Moreover, in contrast 
to the canonical entropy, S{E, N, ..) is everywhere single valued and multiple 
differentiable. There is no need for extensivity, for concavity, for additivity, 
and no need for the thermodynamic limit. This is a great advantage of the 
geometric foundation of equilibrium statistics over the conventional definition 
of the Boltzmann-Gibbs (BG) canonical theory. However, addressing entropy 
to finite eventually small systems we will face a new problem with Zermelo's 
objection against the monotonic rise of entropy, the Second Law. Here the 
Poincare recurrence time might be small and Boltzmann's excuse does not 
work anymore. This is discussed elsewhere (|2l0;0)- 



3 Topological properties of S{E, ■ ■ ■) 

The topology of the Hessian of S{E, ■ ■ ■), the determinant of curvatures of 
S{E,---), determines completely all kinds of phase transitions. This is evi- 
dently so, because e^^^'^~^^'^ is the weight of each energy in the canonical 
partition sum at given T. Consequently, at phase separation this has at least 
two maxima, the two phases. And in between two maxima there must be a 
minimum where the curvature of S{E) is positive. I.e. the positive curvature 
detects phase separation. This is of course true also in the case of several con- 
served control parameters e.g. e = E/L?^n = N/L'^,s = S/L"^, in the case 
of the g = 3 Potts-gas model on a two dimensional lattice of finite size of 
= 50 X 50 lattice points. Here the Hessian is [Aj are the eigenvalues = 
eigen-curvatures] : 
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The whole zoo of phase-transitions can be sharply seen and distinguished 
[fig-(l)]- 
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3. 1 Curvature 

The curvature (Hessian) of S{E,N,- ■ ■) controls the phase transitions see 
ref.(j3). What is the physics behind a positive curvature? For a short-range 
force it is hnked to the interphase surface tension c.f. in fig. 2 and to a negative 
heat capacity. This implys that heat can flow from cold to hot, c.f. fig4. 

3.2 Atomic clusters 

In fig. 3 I show the simulation of a typical fragmentation transition of a system 
of 3000 sodium atoms interacting by realistic (many-body) forces. To compare 
with usual macroscopic conditions, the calculations were done at each energy 
using a volume V{E) such that the microcanonical pressure P = |^/|^ = 
latm. Evidently, the convex region of S{E, P) is the most interesting region. 

3.3 Stars 

Self-gravitation leads to a non-extensive potential energy oc A^^. No thermo- 
dynamic limit exists for E/N and no canonical treatment makes sense. At 
negative total energies often these systems have a negative heat capacity. This 
was for a long time considered as an absurd situation within canonical statisti- 
cal mechanics with its thermodynamic "limit" . However, within our geometric 
theory this is just a simple example of the pseudo-Riemannian topology of the 
microcanonical entropy S{E, N) provided that high densities with their non- 
gravitational physics, like nuclear hydrogen burning, are excluded. We treated 
the various phases of a self-gravitating cloud of particles as function of the 
total energy and angular momentum, c.f. fig.5 and the quoted PRL-paper. 
Clearly these are the most important constraint in astrophysics. In fig. 6 the 
global phase diagram of a rotating, selfgravitating hydrogen cloud is given. 



4 Conclusion 

Entropy has a simple and elementary definition by the area e^^^''^''"'^ of the 
microcanonical ensemble in the 6A^ dim. phase space. Canonical ensembles are 
not equivalent to the micro-ensemble in the most interesting situations: 

(1) At phase-separation ( ^heat engines !), one gets inhomogeneities, and 

a negative heat capacity or some other negative susceptibihty, conse- 
quently: 

(2) Heat can fiow from cold to hot. 

(3) Phase transitions can be localized sharply and unambiguously in small 
classical or quantum systems, there is no need for finite size scaling to 
identify the transition. 
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(4) Also really large self-gravitating systems can be addressed. 



Entropy rises during the approach to equilibrium, AS > 0, also for small 
mixing systems, i.e. the Second Law is valid even if the Poincare recurrence 
time is not astronomically large (0; H; S H) • 

With this geometric foundation thermo-statistics applies not only to exten- 
sive systems but also to non-extensive ones which have no thermodynamic 
limit. More details are discussed in the references, see also my WEB-page 



http://www.hmi.de/people/gross/ 
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Fig. 1. Global phase diagram or contour plot of the curvature determinant 
(Hessian) for a small system, eqn. (2), of the 2-dim Potts-3 lattice gas with 
50 * 50 lattice points, n is the number of particles per lattice point, e is the total 
energy per lattice point. The line (-2,1) to (0,0) is the ground-state energy of the 
lattice-gas as function of n. The most right curve is the locus of configurations 
with completely random spin-orientations (maximum entropy). The whole physics 
of the model plays between these two boundaries. The different regions of positive 
or negative Hessian, i.e. of respectively negative or positive maximum curvature Ai 
correspond to different phases (in the case studied here A2 is everywhere negative): 
At the dark-gray (green in the color version) lines the Hessian is det = 0, this 
is the boundary of the region of phase separation (the triangle APmB) with a 
negative Hessian (Ai > 0, A2 < 0). Here, we have Pseudo-Riemannian geometry and 
a convex entropy. At the light-gray (blue in the color version) lines is a minimum of 
det(e, n) in the direction of the largest curvature (v;^^-Vdet = 0) and additionally 
det ~ 0, these are lines of second order transition. In the triangle APmC is the 
pure ordered (solid) phase with concave entropy and (det > 0,Ai < 0,A2 < 0). 
Above and right of the line CPmB is the pure disordered (gas) phase again with 
concave s(e,n) and (det > 0, Ai < 0, A2 < 0). The crossing Pm of the boundary 
lines is a multi-critical point. Here we have simultaneously Ai = 0,VAi = 0. It is 
also the critical end-point of the region of phase separation. The lighter-gray (red 
in the color version) region around the multi-critical point Pm corresponds to a 
flat, horizontal region of det(e, '^V^Ai ~ and consequently to a somewhat extended 
cylindrical region of s(e,n), see (|j;la); C (red) is the analytically known position of 
the critical point (second order transition) which the ordinary g = 3 Potts model 
(without vacancies) would have in the thermodynamic limit. C is also approached 
by the line of second-order transition (light-gray/blue) in the present small system. 



6 




1 



0.3 



0.5 



0.7 



0/9 



1.1 



1.3 



Fig. 2. The physical origin of convex regions in S{E) for systems with 
interactions of short range. MMMC (0) simulation of the entropy s{e) per atom 
(e in eV per atom) of a system of A'^o = 1000 sodium atoms at an external pressure 
of 1 atm. At the energy e < ei the system is in the pure liquid phase and at e > 63 in 
the pure gas phase, of course with fluctuations which are proportional to the inverse 
negative curvature of s{e). The latent heat per atom is qiat = 63 — 61. Attention: the 
curve s(e) is artifically sheared by subtracting a linear function 25+e*11.5 in order to 
make the convex intruder visible. s(e) is always a steeply monotonic rising function. 
We clearly see the global concave (downwards bending) nature of s{e) and its convex 
intruder. Its depth is the relative entropy loss due to additional correlations by 
internal interfaces. It scales cx A^~^/^. From this one can calculate the surface tension 
per surface atom asurf/Ttr = ^Sgurf * ^o/^surf- For the present example it is 
asurf/Ttr = 5.68 which should be compared to (Jsurf/Ttr = 9.267 for the bulk. 
The bulk value is approximated by CTsurf/Ttr systematically with rising number of 
atoms N [more details c.f.(0;[3)]- The double tangent (Gibbs construction) is the 
concave hull of s(e). Its derivative gives the Maxwell line in the caloric curve T(e) 
at Ttr- In the thermodynamic limit the intruder would disappear and ^(e) would 
approach the double tangent from below. Nevertheless, even there, the probability 
(X e^* of configurations with phase-separations are suppressed by the (infinitesimal 
small) factor g^^^^^ relative to the pure phases and the distribution remains strictly 
himodal in the canonical ensemble. The region ei < e < 63 of phase separation gets 
lost. 
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Fig. 3. Atomic cluster fragmentation in the convex region of S{E). The 
backbending caloric curve T{E) (blue in the color version), left scale. The number 

of fragments Nfr and the number of surface atoms -^e// ~ >2^'i^^ larger 
fragments (green in the color version) can be read off from the right scale. Ttr is 
the "Maxwell" line dividing T{E) into two equal area parts. The inserts above give 
the mass distribution at the various points along the caloric curve T[E). The label 
"4:1.295" means 1.295 quadrimers on average. This gives a detailed insight into 
what happens with rising excitation energy over the transition (convex) region: At 
the beginning (e* ~ 0.442 eV) the liquid sodium drop evaporates 329 single atoms 
and 7.876 dimers and 1.295 quadrimers on average. At energies e >~ leV the drop 
starts to fragment into several small droplets ("intermediate mass fragments") e.g. 
at point 3: 2726 monomers, 80 dimers, ~5 trimers, ~15 quadrimers and a few heavier 
ones up to 10-mers. The evaporation residue disappears. This multifragmentation 
finishes at point 4. It induces the strong backward swing of the caloric curve T{E). 
Above point 4 one has a gas of free monomers and at the beginning a few dimers. 
This transition scenario has a lot similarity with nuclear multifragmentation. It is 
also shown how the total interphase surface area, proportional to A^^^ = ^ • n'^^^ 
with Ni > 2 (Aj the number of atoms in the ith cluster) stays almost constant up 
to point 3 even though the number of fragments (Aj,, = is monotonic rising 
with increasing excitation. 
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Fig. 4. Heat can flow from hot to 

cold. Potts model (here with q = W 
showing a strong first order transi- 
tion), in the region of phase sepa- 
ration. At ei the system is in the 
pure ordered phase, at 63 in the pure 
disordered phase. A httle above ei 
the temperature T = 1/(3 is higher 
than a httle below 63. Combining 
two parts of the system: one at the 
energy ei + 5e and at the temper- 
ature Ti, the other at the energy 
63 — 5e and temperature T3 < Ti. It 
will equilibrize with a rise of its en- 
tropy, under a dropping of Ti (cool- 
ing) and an energy flow (heat) from 
3 1: i.e.: Heat flows from cold to 
hot! Clausius formulation of the Sec- 
ond Law is violated. Evidently, this 
is not any peculiarity of gravitat- 
ing systems! This is a generic situa- 
tion within classical thermodynam- 
ics even for short-range coupling. It 
has nothing to do with long range 
interactions 
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Contour plots and density profiles 
of a rotating, self-gravitating 
A'-body system sliowing the 
formation of a stable double cluster 
(left) and an unstable ring (right) at 
different energies. The 
double-cluster structure illustrates 
the spontaneous breaking of 
rotational symmetry at intemiediate 
energy and high angular 
momentum. 
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Fig. 5. Phases and Phase-Separation in Rotating, Self-Gravitating Sys- 
tems, Physical Review Letters-July 15, 2002, cover-page, by (Votyakov, Hidmi, De 
Martino, Gross) 
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Fig. 6. Microcanonical global phase- diagram of a cloud of self-gravitating 
and rotating system in a spherical container as function of energy and 
angular-momentum. In this mean-field calculation local densities higher than the 

density where nuclear hydrogen burning starts are excluded. Thus, the gravitational 
collapse to only visible stars is followed. Outside the dashed boundaries only some 
singular points were calculated (e.g. see also previous figure). In the mixed phase 
the largest curvature Ai of S{E,L) is positive. Consequently the heat capacity 
or the correspondent susceptibility is negative. This is of course well known in 
astrophysics. However, the new and important point of our finding is that within 
microcanonical thermodynamics this is a generic property of all phase transitions 
of first order, independently of whether there is a short- or a long-range force that 
organizes the system. The importance of using the microcanonical ensemble with 
angular momentum as the second control parameter cannot be overemphasized: The 
kinetic, "'chaotic energy" E^i^ = E — Ep^t — L'^/21 is maximized at large moment 
of inertia /. Similar equilibrium calculations are done in astro-physics with angular 
velocity Q = L/ 1 as control parameter (~ canonical ensemble on the rotating disk), 
whence E^m = E — Epot — / x 17^/2. Then of course, the maximum of the entropy 
is at small moment of inertia and only mono-stars are found. 
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